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INTRODUCTION 


This  final  scientific  report  Jias  been  prepared  by  Dr.  M.W.  Beckstead  of 
Brigham  Young  University's  Department  of  Chemical  Engineering.  It  reviews  the 
progress  from  May  1979  to  December  1982  on  a  program  to  study  various  aspects 
of  solid  propellant  combustion  phenomena.  The  original  program  included  three 
specific  areas  of  study: 

1.  Diffusion  flame  analysis 

2.  Kinetic  mechanisms  of  double  base  propellants  ' 

3. "  Interaction  effects  of  distributed  particle  combustion  in  acoustic 

envi ronments. 


The  diffusion  flame  is  the  controlling  mechanism  in  the  combustion  of 
composite  propellants.  Therefore,  a  thorough  understanding  of  diffusion 
flames  is  prerequisite  to  understanding  solid  propellant  combustion,  and  was 
the  objective  of  one  task  of  the  program.  The  second  task  of  the  program  was 
motivated  by  the  desire  to  better  understand  transient  combustion  such  as 
occurs  in  ODT  and  combustion  instability.  Double  base  propellant  combustion 
is  dominated  by  the  flame  kinetics  which  appear  to  be  more  easily  studied  than 
those  of  composite  propellants.  The  final  task  of  the  program  was  aimed  at 
understanding  the  effect  of  stability  additives  on  combustion  instability. 


During  the  final  year  of  the  contract,  work  on  the  first  two  tasks  was 
discontinued  and  all  work  was  concentrated  on  the  third  task.  Interim  reports 
covering  the  initial  work  have  been  published  in  December  1980  and  February 
1982.  This  report  summarizes  the  work  performed  during  the  final  phase  of  the 
contract  covering  the  time  from  October  1981  through  December  1982.  y» 
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PROGRAM  APPROACH  AND  OBJECTIVES 


A.  REASON  FOR  RESEARCH— HISTORY  AND  BACKGROUND 

Acoustic  suppressants  are  commonly  added  to  low  smoke  and  smokeless 
propellants  to  avoid  the  problem  of  combustion  instability.  However,  the 
basis  for  selecting  the  type,  size  and  concentration  of  suppressant  particle 
for  optimum  performance  is  not  well  founded.  These  suppressants  apparently 
work  by  one  or  more  of  t.iree  mechanisms: 

(1)  energy  loss  due  to  viscous  dissipation  of  drag  forces, 

(2)  modification  of  the  propellant  combustion  response  function,  or 

(3)  energy  interchange  due  to  distributed  combustion. 

Previous  work  by  Diederichsen  [1]  gave  credence  to  the  idea  of  using  a 
modified  Rijke  tube  to  study  the  effects  of  distributed  combustion  of 
particles  on  the  acoustics  in  a  rocket  motor.  His  work  went  so  far  as  to 
roughly  and  very  qualitatively  show  that  some  particle  additives  seemed  to  be 
more  effective  acoustic  damping  agents  than  others.  Our  purpose  was  to  extend 
Diederichsen' s  approach  and  construct  a  burner  in  order  to  get  quantitative  as 
well  as  qualitative  acoustical  data,  and  then  use  these  data  to  interpret  how 
acoustic  suppressants  work. 

B.  EXPERIMENTAL  METHOD— BURNER  AND  ADVANTAGES 

The  principal  objective  of  this  work  was  to  determine  the  effect  that 
distributed  combustion  of  acoustic  suppressants  has  on  combustion 
instability.  The  experimental  basis  for  the  technical  approach  is  a  modified 
Rijke  tube,  a  gas  burner  with  an  attendant  paddle  to  control  (stop  and  start) 
oscillations  (see  Figure  1).  This  burner  is  probably  more  correctly  called  a 
'gauze  burner1  rather  than  a  Rijke  tube  burner,  as  explained  by  Diederichsen, 
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since  the  oscillations  are  apparently  not  driven  by  a  heat  transfer  mechanism 
from  the  gauze  flameholder  to  the  gases  but  by  a  variation  of  the  normal 
burning  rate  and  velocity  coupled  to  an  acoustic  (node  by  the  pressure 
oscillations.  However,  the  acoustics  of  the  system  are  what  would  be  expected 
from  a  Rijke  tube. 

Our  approach  is  to  stu*'y  three  types  of  particles  (very  reactive, 
slightly  reactive,  and  inert)  to  determine  the  effect  of  variable  energy 
release  on  the  acoustic  wave  and  viscous  damping  effects.  This  variable 
energy  release  is  probably  a  function  of  both  the  amount  of  energy  released 
and  the  rate  at  which  it  is  released  (rate  of  reaction  of  the  burning 
particle).  These  effects  can  be  studied  by  comparing  very  reactive  and 
slightly  reactive  particles,  particularly  particles  which  have  very  different 
heats  and  rates  of  reaction.  The  inert  particles  will  allow  us  to  identify 
the  effect  due  to  just  viscous  particle  damping. 

The  following  are  some  of  the  technical  advantages  of  using  this  modified 
Rijke  tube  or  gauze  burner  approach: 

(1)  An  additive  can  be  tested  independent  of  a  solid  propellant  burning 
surface,  allowing  a  separation  of  suppressant  mechani sms  between  the 
modification  of  the  burning  surface  response  and  distributed  combustion 
effects. 

(2)  The  system  can  be  tested  with  and  without  particles,  thus  readily 
evaluating  the  system  effects. 

(?)  Various  types  of  particles  (react,  e  u,  inert)  can  be  tested  under 
identical  conditions  to  be* ter  identify  actual  mechanistic  action. 

(4)  To  perform  tests,  solid  propellants  do  not  need  to  be  mixed,  cured, 
etc.,  thus  reducing  costs  tremendously. 
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(5)  Particle  addition,  frequency,  0/F  ratio,  and  temperature  can  be 
independently  controlled. 

Data  collection  involves  the  fol lowing: 

(1)  pre-combustion  particle  size  distributions 

(2)  gas  flow  rates 

(3)  particle  flow  rates 

(4)  acoustic  pressure  amplitudes  (at  2  or  more  locations) 

(5)  particle  combustion  tracking  (high  speed  movies) 

(6)  post  combustion  particle  collections  and  examination 

(7)  selected  wall  temperatures 

Data  analysis  will  also  involve  the  following: 

(1)  acoustic  growth  rate  measurements 

(2)  acoustic  mode  shape  reconstruction 

(3)  particle  combustion  times 

(4)  viscous  particle  damping  calculations 

C.  THEORETICAL  ANALYSIS 

In  order  to  analyze  the  experimental  data  and  interpret  and  extrapolate 
the  results,  a  theoretical  analysis  of  the  Rijke  burner  is  necessary. 

D.  EXPECTED  RESULTS  AND  FUTURE  WORK 

It  is  expected  that  several  mechanisms  will  be  shown  to  act  together  to 
dampen  acoustic  oscillations,  with  one  mechanism  perhaps  dominating  under 
certain  conditions  and  another  mechanism  dor.  nating  at  other  times.  Future 
work  is  to  complete  the  experimental  data  gathering,  analyze  the  data, 
complete  the  theoretical  work,  and  then  try  to  identify  controlling  mechanisms 
of  acoustic  solid  particle  interaction  with  the  intent  to  apply  this  knowledge 
to  actual  rocket  motors. 


RESULTS 


A.  INTRODUCTION 

The  burner  has  been  constructed,  diagnostic  equipment  has  been  purchased 
and  assembled,  and  checkout  tests  to  characterize  the  burner  limits  have  been 
mostly  completed.  Preliminary  tests  were  run,  and  current  work,  is  examining 
this  preliminary  data  and  making  test  refinements  and  minor  modifications  to 
the  system. 

The  major  accomplishment  has  been  to  extend  Diederichsen' s  approach  and 
modify  the  burner  so  as  to  be  able  to  get  quantitative  information  from  the 
burner  that  should  not  only  enable  us  to  quantify  the  acoustic  effects  of 
different  propellant  additives,  but  to  also  gain  insight  into  the  mechanisms 
by  which  reactive  burning  particles  influence  acoustic  instability. 

8.  BURNER  CONSTRUCTION 

Figure  1  is  a  schematic  of  the  burner  configuration.  It  is  a  gas  burner 
which  will  allow  independent  variations  in  the  0/F  ratio  and  the  flame 
temperature.  Particles  are  added  in  the  gas  inlet  ana  are  ignited  by  the 
flame.  The  flame  holder  is  located  at  least  one-quarter  wavelength  into  the 
burner  to  allow  standing  wave  oscillations  to  develop  (i.e.  a  Rijke  tube). 
High  frequency  pressure  transd'  °rs  are  used  to  measure  the  acoustic  pt esscre 
amplitude,  which  ranges  fvom  0.005  to  0.050  psi ,  and  a  HYCAM  camera  is  used  to 
observe  the  extent  of  particle  burning  by  observation  through  the  quartz 
wails. 

C.  BURNER  TEST  I  NO 

This  burner  has  been  built,  tested,  modified,  and  thoroughly 
characterized  both  at  steady  conditions  as  well  as  at  various  oscillatory 
conditions.  The  frequency  of  the  oscillations,  varied  by  changing  the 
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upstream  and  downstream  burner  lengths,  ranges  from  700  to  1200  Hz.,  and  can 
easily  be  extended  with  slight  modifications  downward  to  300  Hz  or  upward  to 
2000  Hz. 

We  spent  some  time  and  effort  characterizing  the  burner  and  the  following 
results  are  consistent  with  those  of  Diederichsen.  A  gauze  flame,  similar  to 
a  singing  flame  or  a  Rijke  tube  effect,  can  indeed  be  stabilized  on  a  gauze 
flameholder  at  or  near  an  acoustic  velocity  node.  This  nodal  frequency  of 
oscillations  depend  on  both  the  upstream  and  downstream  tube  lengths,  and 
there  is  a  certain  tube  length  below  which  oscillations  cease. 

We  also  found  that  the  frequency  can  be  controlled  very  closely  with  the 
tube  lengths.  The  0/F  ratio  can  be  controlled  with  gas  flowmeters,  and  the 
oscillation  amplitude  can  be  controlled  to  some  extent  by  the  nitrogen  flow, 
which  also  controls  the  flame  temperature. 

D.  PADDLE.  CONSTRUCTION 

The  major  effort  of  the  past  year  was  to  find  a  way  to  control  (stop  and 
start)  the  oscillations,  since  growth  curves  of  oscillations  give  much  more 
valuable  information  than  does  limiting  amplitude  data,  which  is  ail  that  has 
been  previously  measured.  The  method  selected  to  do  this  and  the  major 
modification  to  Diederichsen' s  approach  was  to  put  a  paddle  on  the  top,  open 
end  of  the  burner  (see  Figure  1)  to  suddenly  block  off  the  burnt  gas  flow  and 
turn  an  acoustic  velocity  anti-node  into  a  velocity  node.  This  effectively 
shuts  down  the  oscillations,  and  opening  the  paddle  again  allows  the 
oscillations  to  grow.  The  resultant  growth  curve  .allows  the  measurement  of 
the  linear  cramping  coefficient,  which  is  a  characteristic  of  the  acoustic  wave. 
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E.  PADDLE  TESTING 


This  paddle  has  been  designee,  constructed  and  tested,  and  has  been  shown 
to  damp  out  oscillations,  though  there  remains  some  residu  ’  noise  (from  six 
to  twelve  percent)  of  the  maximum  limiting  amplitude.  Tin's  noise,  however, 
does  not  seem  to  extend  into  the  linear  part  of  the  growth  curve,  and  is 
usually  a  different  frequency.  Therefore,  it  may  not  be  a  serious  problem. 
However,  there  is  a  significant  problem  concerning  the  speed  the  paddle 
moves.  Due  to  the  unexpected  severity  of  the  acoustic  driving  force  (a  ~  500 
sec-*),  the  growth  curve  is  well  developed  before  the  paddle  can  completely 
finish  uncovering  the  tube.  This  means  that  at  present  our  data  are  suspect 
of  being  influenced  more  by  the  paddle  than  by  the  acoustics  of  the  particles 
and  the  system.  Present  efforts  are  attempting  to  further  modify  the  paddle 
to  increase  its  speed  so  that  it  will  completely  uncover  the  tube  in  less  than 
one  acoustic  oscillation  cycle.  Also,  operating  conditions  are  being  varied 
to  find  conditions  that  will  produce  lower  growth  rates. 

F.  PRELIMINARY  RESULTS 

Preliminary  results  with  small  amounts  of  aluminum  particles  (roughly  1- 
2%)  indicate  that  it  should  be  possible  to  get  the  data  desired.  The  addition 
of  particles  does  si gni f icant iy  change  the  linear  growth  rate  as  expected. 
The  limiting  amplitude  of  the  oscillations  is  also  changed.  In  one  case  (see 
Figure  2}  a  reduced  growth  rate  showed  that  the  particles  were  damping  the 
oscillations,  while  the  limiting  amplitude  increased  indicating  that  the 
particles  could  be  driving  rather  than  damping  the  oscillations.  Ihese  mixed 
results  are  no:  really  contradictory ,  in  that  alpha  is  a  measure  of  the  linear 
acoustic  growth  while  the  limiting  amplitude  is  more  dependent  on  a  multitude 
of  non-linear  effects. 


lime  (sec) 


Figure  2.  Comparison  of  typical  acoustic  growth  data  from  the 

modified  Rijke  Burner  The  upper  cur-ve  is  a  gas  only 
firing.  The  loi/er  curve  was  at  the  identical  conditions 
with  -I-2'T  20pm  A1  added. 


As  mentioned  previously,  the  driving  force  is  much  stronger  than 
originally  expected,  giving  rise  to  an  alpha  much  larger  than  is  desirable. 
Current  efforts  based  on  these  preliminary  results  will  be  to  decrease  the 
severity  of  system  driving,  increase  the  paddle  speed,  extend  the  burner's 
range  of  operation,  and  obtain  the  desired  quantitative  results. 

G.  THEORETICAL  ANALYSIS 

The  theoretical  analysis  of  the  Rijke  tube  system  is  underway  concurrent 
with  the  experimental  work.  The  development  of  a  mathematical  model  of  the 
Rijke  tube  began  with  an  extensive  literature  review.  Modern  workers  that 
have  analyzed  the  Rijke  tube  mathematically  include  Neuringer  and  Hudson  (2), 
Carrier  (3),  Jones  (4),  and  Mugridge  (5).  A  completely  general  analysis  of 
the  Rijke  tube  with  particulate  flow  will  require  solution  of  the  conservation 
equations  for  a  multiphase  flow.  However,  as  the  above  authors  performed 
analyses  of  the  Rijke  tube  without  particle  flow,  the  complete  set  of 
conservation  equations  was  not  solved  by  any  of  them.  The  mathematical 
analyses  on  the  Rijke  tube  performed  by  Neuringer  and  Hudson  (2)  and  Carrier 
(3)  are  typical  of  all  of  the  analyses  performed,  and  are  very  similar. 

The  main  objective  of  the  mathematical  modeling  efforts  of  the  above 
mentioned  authors  was  calculation  of  the  characteristic  frequencies  of  the 
Rijke  tube  and  the  growth  constants  associated  with  each  frequency. 

The  following  assumption^  were  made  by  most  authors: 

1)  All  properties  (pressure,  temperature,  density,  velocity,  thermal 
conductivity,  heat  capacity,  and  ratio  of  specific  heats)  were  taken  as 
constant  in  each  tube  section,  with  a  step  change  in  properties  occurring 
across  the  flame  or  gauze. 

2)  One  dimensional  flow  was  assumed. 

J)  The  effects  of  friction  at  the  tube  wall  were  neglected. 
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4)  Small  amplitude  oscillations  were  assumed,  thus  allowing 
linearization  of  the  conservation  equations. 

5)  The  acoustic  oscillations  were  assumed  to  be  sinusoidal  with  respect 
to  time. 

6)  The  flame  or  heater  gauze  was  taken  as  thin  enough  to  be  neglected, 
and  in  most  cases  was  treated  as  a  simple  discontinuity  between  the  cold  and 
hot  tube  sections. 

7)  Most  analyses  neglected  the  effects  of  gravity  and  other  body  forces. 

8)  The  flame  speed  was  assumed  to  be  negligible  compered  to  the  speed  of 
sound. 

The  primary  difference  between  the  models  was  the  manner  in  which  the 

fluctuating  heat  transfer  from  the  gauze  or  flame  was  represented.  As  the 

continuity  conditions  of  the  gauze  or  flame  were  derived  from  the  mathematical 
equations  chosen  to  describe  the  heat  source,  the  values  obtained  for  the 

characteri stic  frequencies  and  associated  growth  constants  were  effected  by 
the  flame  representation. 

As  a  result  of  the  similarity  between  Carrier's  approach  and  that  of 

Neuringer  and  Hudson,  the  final  eigenvalue  equation  obtained  by  Carrier  was 
quite  similar  to  that  obtained  by  Neuringer  and  Hudson.  However,  there  were 
some  differences  between  the  two  analyses. 

The  most  significant  difference  was  that  Carrier  did  not  apply 
assumptions  2)  and  3),  but  used  two-dimensional  cylindrical  coordinates  in  his 
analysis.  Thus,  the  linearized  partial  differential  equations  obtained  were 
functions  of  the  radial  coordinate  as  well  as  the  axial  coordinate  and  time. 
The  radial  portion  of  the  solution  accounted  for  the  effects  of  viscous 
dissipation  at  the  tube  walls.  The  effect  of  this  was  a  small  but  significant 
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change  in  the  values  of  the  complex  eigenvalues  accounting  for  the  additional 
damping  of  oscillations  due  to  viscous  drag  at  the  tube  walls. 

Second,  Carrier  modeled  the  gauze  as  a  series  of  parallel  heated 
strips.  Using  approximate  boundary  layer  equations,  he  obtained  an  equation 
for  predicting  the  variable  heat  transfer  rate  from  the  heated  strips  to  an 
oscillating  fluid.  His  attempt  to  model  the  heat  release  rate  front  the  gauze 
represented  an  improvement,  though  not  a  significant  difference  between  the 
two  analyses. 

Finally,  Carrier  attempted  to  formulate  the  open  end  boundary  conditions 
more  accurately  than  Neuringer  and  Hudson  by  drawing  upon  acoustic  theory.  Of 
the  three  differences  between  the  analyses,  the  third  is  probably  the  least 
significant. 

Although  no  articles  were  found  in  which  operation  of  the  Rijke  tube  with 
particle-laden  flow  was  discussed  specifically,  some  authors  have  discussed 
the  effect  of  particles  on  oscillating  flows. 

Dobbins  and  Temkin  (6)  and  Culick  (7)  both  analysed  mathematically  the 
effects  of  particles  carried  in  an  oscillating  stream  of  gas.  Dobbins  and 
Temkin  obtained  an  expression  giving  both  the  attenuation  of  the  wave  and  the 
dispersion  (change  in  the  speed  of  sound)  due  to  the  particles.  Culick  (7) 
formulated  a  generalized  equation  for  calculating  the  growth  constants  and 
frequencies  of  longitudinal  oscillations  in  a  chamber  of  variable  cross- 
section  with  particulate  flow  and  mass  addition  from  the  perimeter.  He  also 
demonstrated  the  use  of  his  general  equations  by  investigating  the  dynamics  of 
the  Rijke  tube. 

Expanding  upon  his  earlier  work,  Culick  (8,9)  examined  the  influence  of 
residual  combustion  by  the  particles  in  a  particle-laden  flow  on  the  frequency 
and  attenuation.  In  these  works,  Culick  accounted  for  the  effects  of 
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combustion  due  to  transformation  of  the  particles  from  one  composition  to 
another,  but  he  did  not  account  for  the  effect  of  the  i nteractu^  between  the 
actual  flame  surrounding  the  burning  particle  and  the  oscillating  gas. 

Smoot  and  Pratt  (10)  discussed  the  formulation  of  the  conservation 
equations  for  multicomponent  flow.  Their  specific  application  was  the 
combustion  of  pulverized  coal  dust  in  air,  but  the  equations  formulated  were 
designed  to  handle  any  reacting  or  non-reacting  particle-laden  flow.  These 
general  equations  will  be  the  start' ng  point  for  the  work  to  be  undertaken. 
The  analysis  to  be  performed  will  include  the  damping  of  oscillations  by 
particles  of  arbitrary  size  distribution  in  a  Rijke  tube  having  a  varying 
temperature  distribution  in  the  hot  gas  section.  This  modeling  effort  will 
apply  Culick's  work  to  the  equations  of  Smoot  and  Pratt.  Integration  of  the 
resulting  differential  equations  will  be  accomplished  by  combination  of 
analytical  and  numerical  integration. 

The  objective  of  the  modeling  effort  is  to  predict  the  frequencies  and 
growth  constants  of  oscillations  producted  by  a  Rijke  tube,  both  with  and 
without  particle-laden  flow.  An  attempt  will  be  made  to  account  for  the 
effects  of  the  combustion  of  the  particles. 

Detailed  modeling  of  the  response  of  the  flame  to  pressure  fluctuations 
will  be  avoided,  as  this  is  beyond  the  scope  of  the  present  study. 
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